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ABSTRACT 

Context. Interstellar dust is a key element in our understanding of the interstellar medium and star formation. The manner in which 
dust populations evolve with the excitation and the physical conditions is a first step in the comprehension of the evolution of inter- 
stellar dust. 

Aims. Within the framework of the Evolution of interstellar dust Herschel key program, we have acquired PACS and SPIRE spec- 
trophotometric observations of various photodissociation regions, to characterise this evolution. The aim of this paper is to trace the 
evolution of dust grains in the Orion Bar photodissociation region. 

Methods. We use Herschel/FACS (70 and 160 fim) and SPIRE (250, 350 and 500 fim) together with Spitzer/IRAC observations to 
map the spatial distribution of the dust populations across the Bar. Brightness profiles are modelled using the DustEM model coupled 
with a radiative transfer code. 

Results. Thanks to Herschel, we are able to probe finely the dust emission of the densest parts of the Orion Bar with a resolution 
from 5.6'' to 35. T'. These new observations allow us to infer the temperature of the biggest grains at diff'erent positions in the Bar, 
which reveals a gradient from ~ 80 K to 40 K coupled with an increase of the spectral emissivity index from the ionization front to 
the densest regions. 

Combining Spitzer/IRAC observations, which are sensitive to the dust emission from the surface, with Herschel maps, we have been 
able to measure the Orion Bar emission from 3.6 to 500 jim. We find a stratification in the diff'erent dust components which can be re- 
produced quantitatively by a simple radiative transfer model without dust evolution (diff'use ISM abundances and optical properties). 
However including dust evolution is needed to explain the brightness in each band. PAH abundance variations, or a combination 
of PAH abundance variations with an emissivity enhancement of the biggest grains due to coagulation give good results. Another 
hypothesis is to consider a diff'erent length along the line of sight of the Bar at the ionization front than in the densest parts. 

Key words. Infrared: ISM - ISM: Dust - ISM: photon-dominated region - ISM: individual (Orion Bar) - Radiative transfer 



1. Introduction 

Interstellar grains constantly interact with their surrounding 
gaseous environment. They also strongly interact with radiation, 
absorbing the UV/visible photons from stars and re-emitting 
the energy in the infrared range. Their emission wavelength 
range depends on both grain size and composition. To study 
the dust properties and their evolution, many models have been 
developed in the last two decades in order to reproduce the dust 
emission spectrum and the extinction curve (e.g., Desert et al. 
1990; Zubko et al. 2004; Draine & Li 2007; Compiegne et al. 
2011). Most of these use three dust components to reproduce 
the observational constraints: poly cyclic aromatic hydrocarbons 
(PAHs), "very small grains" (VSGs), and big grains (BGs). 
The two first are made of carbonaceous particles. Due to their 
sizes (~ [0.4 - 10] nm), they are stochastically heated and emit 
most of their energy below 60 jim. The BGs (up to ~ 200 nm) 
are in thermal equilibrium with the radiation field and contain 
most of the dust mass. During their life, dust grains undergo 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



modification of their size, structure and chemical composition. 
To understand this evolution, we need to trace the characteristics 
of grains in relation to changes of the physical and dynamical 
properties. 

Photo-dissociation regions (PDRs) are the surface layers where 
the radiation field is able to dissociate H2 molecules but 
cannot ionize hydrogen atoms. Thus, they are found at the 
interface of HII regions and molecular clouds. In such zones, 
the radiation field coming from nearby OB stars regulates the 
physical and chemical evolution of the gas. That is why PDRs 
are the favourite laboratories to study evolution of the dust 
and gas constituents with the local excitation and the physical 
conditions. 

The Orion Bar has been extensively observed, leading to a large 
dataset in diflTerent tracers of the gas component (e.g., White 
& Sandell 1995; van der Werf et al. 1996; Walmsley et al. 
2000; Goicoechea et al. 2011). Located between the Trapezium 
cluster and the Orion Molecular Cloud, the Bar is part of the 
Orion nebula lying at 414 ± 7 pc from the Earth (Menten et al. 
2007).The incident radiation field is dominated by the 06 star 
Ori C and has been estimated to G = [1 - 4] x 10"^ Go at 
the ionization front (Marconi et al. 1998; Tielens & Hollenbach 
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Fig. 1. Orion bar maps observed by IRAC, PACS and SPIRE instruments. The red circle bottom right stands for the EWHM for each 
channel, and the black star shows the location of the illuminating source. The saturated pixels are indicated in black. 



1985) where Go is the integrated intensity of the standard 
interstellar radiation field (ISRE) given by Habing (1968). 
Observational studies of the bright Bar have been accompanied 
by several modelling studies to determine its physical structure, 
excitation, and gas emission. Previous studies revealed the pres- 
ence of diff'erent layers within the PDR, with a stratified structure 
typical of an edge-on geometry (e.g., Tielens et al. 1993; Tauber 
et al. 1994; Lis & Schilke 2003). HCO+ and HCN / = 1 - 
emission line maps from Young Owl et al. (2000) showed the 
presence of dense clumps (jih = 3x10^ cm"^) embedded in an 
interclump medium with a density nn ~ 5 x 10^ cm"^, which 
emits most of the intensity. 

Within the framework of the "Evolution of interstellar dust" 
Herschel key program (Abergel et al. 2010), the Orion Bar 
has been observed by the ESA Herschel Space Observatory 
(Pilbratt et al. 2010). Its unprecedented spatial resolution to- 
gether with its spectral coverage at the far-Infrared wavelengths 
allow us to probe the dust evolution within PDRs. The large sen- 
sitivity range and the spectrophotometric capabilities of PACS 
(Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) detectors 
give a unique view of the Orion Bar at the far-Infrared wave- 
lengths. Several papers by our group deal with Herschel obser- 
vations of the Orion Bar. Bernard-Salas et al. (2012) focus on 
PACS spectroscopic data whereas Habart et al. (in prep.) explore 
SPIRE/FTS observations. In this paper, we present the photo- 
metric maps from PACS and SPIRE to study the emission of 
dust at thermal equilibrium with the radiation. We seek to ad- 
dress the following issues: what are the processes involved in 
the dust evolution in the Orion Bar; are there abundance varia- 
tions from the diff'use ISM; and are there any diff'erences in the 
optical properties of the grains in the Bar compared to grains in 
the diff'use ISM? 

Section 2 describes the photometric data from SPIRE and 
PACS, explaining how they have been reduced. In Section 3, the 



morphology of the Orion Bar in Spitzer and Herschel bands is 
presented. These new observations allow us to study both the 
evolution of the BG spectral energy distribution across the Orion 
Bar (Sect. 4), but also to probe dust emission by mapping its spa- 
tial distribution (Sect. 5). In section 6, we model the dust emis- 
sion in the Orion Bar using the DustEM model (Compiegne et al. 
201 1) coupled with a radiative transfer code, which is compared 
to our data. The results are discussed in section 7. 



2. Observations and data processing 

The Orion Bar was mapped during the Herschel Science 
Demonstration Phase (SDP) on February 23 and March 1 1 2010 
using the PACS (70 and 160 yum) and SPIRE photometers (250, 
350 and 500 fim). For PACS, two concatenated and perpendicu- 
lar 12'xl2' scan maps (corresponding to two obslds) were ob- 
served using the medium scan speed (20"/s), a scan length of 
12', a cross-scan step of 50", and a 15 scan legs (total observ- 
ing time 1726 s). For SPIRE, a single map was taken without a 
repetition (total observing time of 193 s). 



2.1. Data processing 

The SPIRE maps reported in this paper are the Level 2 naive 
maps delivered by the Herschel Space Center (HIPE version 
7.0.1991), with standard corrections for instrumental eff'ects 
and glitches. Striping induced by off'sets in the flux calibration 
from one detector to another was removed using the Scan 
Map Destriper module included in the HIPE environment. The 
overall absolute flux accuracy is dominated by the calibration 
uncertainty and is conservatively estimated as ±7%^. 



^ from the SPIRE observers' manual, available on the webpage: 
http : //Herschel . esac . esa . int/ 
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The PACS data were processed with HIPE (version 6.0.1196) 
with special care for deglitching and the removal of the 1/f 
noise component. Eor each filter, processing from Level to 
Level 1 is performed on each obsid and yielded one Level 1 
frame by obsld. A second level glitch mask was computed from 
preliminary maps built from the frames after high-pass median 
filtering. Object masks were also attached to the Level 1 frames. 
These were computed from the preliminary maps (using the 
second level glitch mask) to flag out the brightest parts of the 
sky on the high-pass filtering step performed from Level 1 to 
Level 2 (see below). 

The processing from Level 1 to Level 2 combines the two 
Level 1 frames in order to obtain one final map for each filter. 
The 1/f noise component was removed with high-pass median 
filtering, using the glitch mask and the object mask. In order 
to avoid any filtering of extended emission, we have taken a 
relatively broad window size, corresponding to 5 times the scan 
length. Einally we performed a simple coaddition of the data. 
The overall absolute flux accuracy is estimated to be ±20% 
(Ali 2011) for both the blue and red bands corresponding to the 
conservative value in the case of extended sources. 
The maps from our Herschel programme are shown in Eigure 1 
with the 3.6 //m IRAC observation (from the Spitzer Data 
Archive). The zero point of Herschel maps is unknown there- 
fore, a background emission centered around the position 
= 5h 35min 26.7s and dj^^^^^ = -5°26'4.7'' and corre- 
sponding to emission unrelated to the Bar was subtracted from 
the data. This point is not critical for our analysis given the 
high brightness of the Orion Bar compared to the background 
emission. The five processed maps are publicly available at 
http://idoc-herschel.ias.u-psud.fr/sitools/client_user/. 

2.2. Beam effects 

To make a coherent spatial study and extract spectral energy dis- 
tributions (SEDs), all the maps need to be brought in the same 
beam, in our case the 500 jim beam. The simplest way to do that 
is to assume that Herschel beams are Gaussian. Therefore each 
map is usually convolved by a Gaussian which full width at half 
maximum (EWHM) is given by: 



EWHM = yEWHM^ - EWHM^, (1) 

where EWHMi is the EWHM of the considered beam and 
EWHM2 is the EWHM of the larger beam of the dataset. 
Nevertheless, simulations and measurements of Herschel point 
spread functions (PSEs) show that the real beams are signifi- 
cantly diff'erent from Gaussian since they present asymmetries 
and secondary lobes (Sibthorpe et al. 2011). This point is partic- 
ularly critical since ignoring secondary lobes misestimates the 
fluxes at the edges of bright structures as illustrated in Eigure 2. 
Thus, it is necessary to take into account the real shape of the 
beam in our study and we use the theoretical PSEs delivered by 
the Herschel Science Centre. The transition functions needed to 
convolve the maps to a common beam size does not have an an- 
alytical form, and must be numerically calculated. 

Computing a transition function k between two PSEs, 
PSE2 = PSEi ^, is a deconvolution problem (Gonzalez & 
Woods 1992). We choose a regularized least- square method 
(Tikhonov 1963) to compute a transition PSE k from a given 
channel of PSEi to PSE2: 

k = SLVgrninJik) = IIPSE2 - PSEi ^ k\\l -h /i \mk)\\l , (2) 



% 




Fig. 2. (/250/im * h(x, FWHM) - hso^m * ^250->50o): relative difl'er- 
ences between the 250yum observation brought to the 500 yum 
resolution using a 30.4'' EWHM Gaussian and the transition 
PSE obtained with our method; black contours show the 250 yum 
emission (levels: 10500, 14000, 17500 MJy.sr'^). 
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Fig. 3. Relative error in the reconstruction of the SPIRE 
500 yum from the 250yum one using our inversion method: 
(PSE500 - PSE250 ^ k25o->50o)/max(PSE5oo). 



where Q is a smoothness function taken as a derivative in our 
case, and jd is the regularization parameter which can be fixed 
according to the noise level present in both PSEs. A transition 
PSE corresponds to a couple of channels for a given observa- 
tion because the orientation angles of the instruments are re- 
quired in the computation. Recently, Aniano et al. (2011) have 
proposed common-resolution convolution kernels based on an- 
other method. In their approach, the low pass filtering, needed 
to avoid problems in the high frequency component, is made by 
a filter which shape is fixed whatever the considered beams. In 
our method the filtering is optimized by adapting the filter ac- 
cording to the shape of the considered PSEs, thus the transition 
function computed is optimal in the least-squares sense. Thanks 
to regularization process, this method is able to give a very good 
reconstruction of respecting the balance between spectral in- 
formation available in the PSEs and the noise levels (Eig. 3). We 



3 



H. Arab et al.: Evolution of dust in the Orion Bar with Herschel 



IRAC 5,6 um 




145.00 



0.00 



0.00 



0" 

23- 30" 

0" 

24- 30" 

0" 

25- 30" 

0" 

26- 30" 

0" 

-5° 27- 30" 



SPIRE 250 um 



'9E+04 



9.5E+03 



0.00 




+04 



3.2E+04 



SPIRE 500 um If^r/9^+03 



550.00 



0.00 



20^ 15^ 
J2000 R. A. 



10^ 



Fig. 4. Observed maps convolved to the 500 yum beam (red circle; FWHM=35.1") using the inversion method described section 2.2. 
The three colored squares on the 500 fim observation stand for the positions of the three SEDs shown Eigure 5. Black contours show 
the 250)um emission convolved to the 500 //m beam (red circle) with levels at 10000, 12500 and 15000 MJy.sr"^ 



present the Herschel observations convolved to the 500 jim res- 
olution and gridding in Eigure 4. 



3. Morphology of the Orion Bar with Herschel 

We have mapped the emission of the BGs, which trace the 
matter up to densest regions, of the Orion Bar. Indeed, with 
Spitzer/IRAC, we only detect the emission of the smallest dust 
particles which is strongly dependent on the radiation field 
since they are stochastically heated. Therefore, we observe the 
emission coming from the surface of dense regions. This is 
illustrated by the first panel of Eigure 1. The 3.6 jjm map from 
Spitzer/IRAC exhibits numerous narrow bright filaments due to 
the illuminated edge of the Bar. These bright filaments contain 
small scale structures, revealing the complexity of the edge of 
the Bar. 



Spitzer/IRAC data, we are able to study the spectrum of the BG 
emission and to make a coherent spatial study of the diff'erent 
dust populations in the Orion Bar. 



Table 1. Contribution of VSGs to the overall intensity (in %) 
according to the DustEM modeF in the five Herschel bands for 
the spectra shown in Eig.5. 



Offset from the 500 fim peak posi- 
tion (5 h 35 min 16 s;-5°23'23'') 


10/um 


\60/um 


250/um 


350/dm 


SOOjim 


(+2S'\-14'') (top panel of Fig 5) 


4.1 


3.5 


3.6 


4.0 


5.4 


(0",0") (middle panel of Fig 5) 


4.0 


2.4 


2.3 


2.6 


2.7 


(-28",+28") (low panel of Fig 5) 


4.2 


2.0 


1.9 


2.1 


2.2 



Notes. We use the same dust populations as in Compiegne et al. 
(2011) and an excitation corresponding to 6^ Ori C. 



In Herschel bands, it is striking to note that the global 
structure of the Bar appears the same (Eig. 1), especially if we 
look at the maps brought to the same resolution (Eig. 4). Even 
if the ridge appears broader because of resolution eff'ects, the 
shape is the same and we recognize the sub-structure along the 
Bar seen with IRAC. However, as we are sensitive to BGs and 
denser matter, we see several diff'erences in the sub- structure. 
Eor example, the South- West part of the Bar ends with a bright 
cloud invisible at shorter wavelengths, indicating the presence of 
cold material. Besides, the sub- structures detected along the Bar 
have not the same relative brightness at diff'erent wavelengths. 
The brightest one at 70 jim is the faintest at 250 yum, revealing 
a strong difference of temperature. Finally, we observe with 
Herschel an increase of the distance between the observed Bar 
and the illuminating star with increasing wavelength. 

With these new Herschel observations, combined with 



4. BG spectrum 

The PACS and SPIRE data allow us to derive the BG component 
spectrum at each position on the Bar. After convolving all maps 
to the 500 fim beam, we fit a modified blackbody to the spectra 
extracted according to: 

Iy=Ty^{^^j B,(n (3) 

where ly is the specific intensity, Ty^ is the dust optical depth at 
frequency yo,y5 is the spectral emissivity index. By is the Planck's 
function, and T is the dust temperature. The fits are performed 
using the MPFIT IDL function (Markwardt 2009), which re- 
lies on the Levenberg-Marquardt algorithm. The free parameters 
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Typ, p and T are computed for each pixel of the maps. The five 
bands are included in the fits since BGs dominate in all bands. 
This has been verified with DustEM, which is able to quantify 
the contribution of non-equilibrium emission at a given wave- 
length (Tab. 1). 

Errors entered in the fits are conservative and correspond to 
the calibration uncertainties (20% for PACS and 7% for SPIRE) 
since statistical noise is negligible compared to calibration er- 
rors. The fitted spectra are integrated over the instrument filters 
so as to take the colour correction into account. We can then 
study the behaviour of the EG spectrum at diff'erent positions 
in the cloud. Figure 5 shows the evolution of the BG spectrum 
while crossing the PDR. The exact positions are shown by the 
colored squares in Figure 4. In front of the bar (top panel of 
Fig. 5), the fitted temperature is 81 ± 25 K whereas as we en- 
ter into the bright ridge, it decreases sharply to 47.5+ 7 K (mid- 
dle panel of Fig. 5). Behind the Bar (low panel of Fig. 5), BGs 
cool down to 39+4 K. The p spectral emissivity index follows 
an inverse behaviour, increasing from 1.10 to 2.08. This trend is 
real but might be amplified by the mixture of grains of diff'erent 
temperature along the line of sight. Indeed, mixing can broaden 
and flatten significantly the modified blackbody shape. This ef- 
fect is more important in front of the Bar, where the geometry is 
not edge-on (purple spectrum). Thanks to the high resolution of 
Herschel observations, we are able to probe T and p variations 
on small spatial scales (~35'0, in a single and transition object 
for the first time. A full study of T and p in the Orion Bar region 
will be presented in detail and compared to previous results in a 
second paper (Arab et al. in prep). 



5. Brightness profiles 

Using the Herschel data coupled with Spitzer/IRAC observa- 
tions, we have been able to study the spatial variations of the dif- 
ferent dust population. The BG emission is probed by Herschel 
bands whereas the 3.6 yum emission from IRAC is used as a 
proxy of the PAHs. The 8 jim channel, usually taken as a tracer of 
aromatic emission, is aff'ected by saturation. Besides, Herschel 
resolution off'ers the opportunity to resolve the BG emission 
within the Orion Bar. 

We extracted brightness profiles obtained along a cut going 
from the exciting star through the PDR, whose position is vis- 
ible in Figure 7. The Orion Bar appears as a prominent peak 
in these cross-cuts lying at ~ T from 6^ Ori C (Fig. 6). Sub- 
structure in this profiles appears at 3.6 and 70 yum showing that 
the Bar is more complicated than a simple ridge. A shift be- 
tween the peak positions is also clearly visible. This shift is real 
as can be seen on the comparison of the 250 yum map superim- 
posed with contours of the 3.6 yum map after convolution and 
reprojection as described above (Fig. 7). The 250yum emission 
clearly rises further from the exciting star than the 3.6 yum one. 
This can be generalized: the shorter the wavelength is, the closer 
to the exciting source is the emission. As we probe bigger grains 
at larger wavelengths, we see that the BG emission is located 
further the ionization front than the emission from the smaller 
grains. Moreover, we show in Figure 7 the C^^O emission ob- 
served by the SPIRE-FTS (Habart et al., in prep.), which traces 
the dense matter and a good correlation is observed with the 
250 yum broadband emission, confirming that the FIR dust emis- 
sion is a tracer of the densest part of the PDR (see also Buckle 
et al. 2011). 

The ionization front is located at 1 1 from 6^ Ori C (Pellegrini 
et al. (2009) from [SII] intensity profile; Bernard-Salas et al. 



10* 



0) 

c 



10=^ 







I In front of the Bor 


T=70.6±10.5 K : 




p= 1 .ZZ±(JJ4 . 


r / 






1 ■ 



10 



100 

Wovelength (/zm) 



1000 



10* 



(n 
c 

(U 



10^ 





■ 1 ' : 


: Inside the Bor 


T=48,8± 4.0 K : 




j8=1.62±0J5 : 











10 



100 

Wovelength (/xm) 



1000 



10* 



tn 
c 

(D 



10^ 





■1-1 1 1 1 — 1 — I I I I ■ 


: Behind the Bor 


T=37J ± 2.5 K : 




j8 = 2.17±0.15 : 




xl =1^3^ 




■ 1 ■ \ ■ 



10 



100 

Wovelength (/znn) 



1000 



Fig. 5. Modified blackbodies fits at diff'erent positions across the 
Bar. Reduced is indicated on each fit. The position of each 
SED is indicated in the 500 yum panel of Figure 4. 



(2012) from [Nil] PACS observations). Adopting the same cut 
as Bernard-Salas et al. (2012), we find from the 3.6 jim emission 
that PAHs peak at 1 16 +2'' from the star, just behind the ioniza- 
tion front. We lack good VSG tracers since 24 yum Spitzer/MIFS 
data of this region are unusable because of saturation and the 70 
yum PACS emission is dominated by BGs. The peak of the BG 
emission is around 150 far from the exciting star. The fact that 
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Fig. 6. Normalized brightness profiles across the Orion Bar at 
native resolution (black solid lines) and convolved to the 500 yum 
(grey dashed lines). The resolution of the solid line profiles cor- 
responds to the nominal resolution of each instrument: 1.88'' for 
IRAC, 5.6'' for PACS 70 yum, 11.3." for PACS 160 yum, 17.6" 
for SPIRE 250 yum, 23.9" for SPIRE 350 yum, 35.1" for SPIRE 
500 yum. 
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Fig. 7. Herschel 250 yum observation (color image) and IRAC 
3.6 yum map (red contours) in the SPIRE 250 fim beam (red 
circle; FWHM=17.6"). Black contours show the distribution 
of C^^O (8-7) observed with the SPIRE-FTS (yellow circle; 
FWHM=23.9") from Habart et al. (in prep.). The white star il- 
lustrates the position of 6^ Ori C and the white line is the cut of 
Figure 6. 



the emission of diff'erent species rises successively when we en- 
ter the PDR is a result of the nearly edge-on orientation of the 
Orion Bar (see also Dicker et al. 2009). 



in the Orion Bar with Herschel 
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Fig. 8. Density profile with a power-law rise and a symmetric 
fall. The maximum density is 1.5x10^ cm"^ and is reached at 
5.25x10-2 pc far from the PDR edge. 

6. Modelling 

The Orion Bar is a very dense PDR which is illuminated by a 
strong radiation field, therefore the dust size distribution and op- 
tical properties should show some diff'erences to those seen to- 
ward diff'use regions. Our goal is to extract evidence of this evo- 
lution, and to characterise and to quantify it. However, radiative 
transfer eff'ects have a crucial influence on the dust emission, 
which is why modelling is needed to recover the properties of 
grains in the Orion Bar. 

6.1. Model description 

We model the brightness profiles using the DustEM code for 
the dust emission coupled with a plane parallel radiative transfer 
code. 

DustEM (Compiegne et al. 2011) is a numerical tool which cal- 
culates the emission of an ensemble of dust populations of vari- 
ous properties under the eff'ect of a given radiation field. The dust 
model we use includes four diff'erent populations: PAHs, small 
amorphous carbons (SamC), large amorphous carbons (LamC), 
and astronomical silicates (aSil). The dust abundances and emis- 
sivities used in the model are indicated in Tab. 2. The reader is 
referred to Fig. 1 of Compiegne et al. (201 1) to get a detailed de- 
scription of the size distribution and to Fig. 4 and 5 in the same 
reference to examine the computed extinction curve. In PDRs, 
evolution of the physical conditions is dominated by the changes 
in the radiation field. Therefore the radiative transfer eff'ects are 
critical for our study and need to be properly quantified. We have 
used the model described in Compiegne et al. (2008) to take 
them into account. In this model, the PDR is represented by a 
semi-infinite plane-parallel slab whose density is defined by the 
user. The cloud is illuminated by the Mathis ISRF (Mathis et al. 
1983), plus a blackbody radiation field derived from stellar pa- 
rameters of 0^ Ori C (radius and distance to the PDR). Then the 
code computes the flux F(z, A)\ 

F(z, A) = Ftiz, A) + Fbiz, A) + Fir(z, A), (4) 

where z is the depth in the PDR, A is the wavelength, Ft(z, A) 
is the transmitted flux, Fiy(z,A) is the backscattered one, and 
Fir(z,A) denotes the emission coming from low depth dust. 
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which absorbs, re-emits and thus heats the grains located deeper 
in the PDR. The calculation was performed assuming an asym- 
metry parameter for the scattering phase function g of 0.6 which 
corresponds to the generally agreed value for the diffuse ISM 
(Witt et al. 1997). Moreover, the model only considers a single 
backscattering, in other words, when a photon is backscattered, 
it cannot go back toward the densest parts of the PDR. This ap- 
proximation is justified as the g value indicates a strong forward 
scattering (88% forward against 12% back). Nevertheless, the 
main assumption of this model is the plane-parallel approxima- 
tion. It greatly simplifies the computations since it allows a 1-D 
treatment. 

The radiative code requires a density profile hh^z) for the 
PDR. For convenience, we use a power-law rise starting at zo 
and defined as: 

nniz) = noX |^ j , (5) 

where zo is the depth at which the density reaches the maximum 
value no and a is the power-law index. The density keeps 
constant at the maximum value between zo and zi and then 
decreases symmetrically. Indeed, models with a rise followed 
by a constant density do not reproduce the observed decrease of 
intensity after the peak at larger wavelengths. A lower density is 
required and for convenience we adopt a symmetric fall off. The 
density profile is showed in Figure 8. 

The four parameters, no, a, zo and zi are free to vary in the 
fitting, contrary to the radiation field parameters (the radius of 
the exciting star, its temperature, and its distance from the cloud) 
which are fixed to give an illumination equal to 2.6 x 10"^ Go 
in Mathis units, following Marconi et al. (1998). Moreover, as 
the DustEM code computes the emission per hydrogen atom it 
is thus necessary to make an assumption on the column density 
to compare our model outputs to the data. It is made through 
the IpDR parameter, corresponding to the PDR length along 
the line of sight. It is related to the column density Nntotiz) by 
Nntotiz) = nniz) Ipdr, where n(l) is the density along the line 
of sight. This is a fifth parameter which is set to reproduce the 
emission given a density profile. Ipdr has to be compatible with 
previous measurement of the column density in the Orion bar. 
2x10^^ cm"^ is an order of magnitude given by Johnstone et al. 
(2003) from SCUBA/JCMT observations at 850 yum. 

To take into account the self-absorption by dust, which can 
be important in the NIR, the emerging intensity is reduced by: 

Hpdr p-T{A,l) 

dl, (6) 

Jo ^PDR 

where / is the distance along the line of sight (/ = at the PDR 
edge) and t(/1, /), the optical depth, computed by assuming no 
density gradient along the line of sight and using the extinc- 
tion curve computed by DustEM. This factor is only significant 
at 3.6 yum since the optical depth is actually very low at longer 
wavelengths (Tmax < 0.4 at 70 yum). Finally, the model output 
is convoluted by the appropriate PSF and integrated over each 
considered filter to be compared to the data. 

6.2. Fitting 

Several different density profiles have been tested to explore 
the parameter space and the influence of each parameter on the 
dust emission profiles. We allow no to vary between 10^ and 
10^ cm"^ which is the canonical density range for the Orion Bar 



(Hogerheijde et al. 1995; Tauber et al. 1994). For each no value, 
several power-law indexes a are considered between 1 (a linear 
rise of density) and 12 (a very steep profile). As for zo and zi, 
they are well constrained by the maximum intensity position in 
each band and require to be fitted a posteriori given the values 
of no and a. 

The maximum density no affects the overall brightness. The 
densest is the PDR, the highest is the brightness in Herschel 
bands. For a given no and at a given wavelength, the a power- 
law index defines where the dust emission peaks. Indeed, the 
brightness peak position depends on a competition between the 
incident flux and the density of absorbing particles. There is thus 
a degeneracy between no and a and to keep a peak position 
constant with an increasing density, a higher power-law index is 
needed. We also notice that the couple (^o. changes the shape 
of the profiles. On the other hand, the length of the PDR along 
the line of sight is a multiplicative factor and affects each band 
in the same way. 

Concerning the dust model, we consider the abundances and 
optical properties of dust in the diffuse ISM. The best adjust- 
ment of the data has been obtained with the following param- 
eters: no = 1.5 X 10^ cm"^, a = 2, zo = 5.25 x 10"^ pc, 
Zi = 6.95 X 10"^ pc and /pdr = 0.45 pc. Using these param- 
eters we obtain Nh = no x I pdr = 2 x 10^^ cm"^ which is of 
the same order of magnitude as measured by Johnstone et al. 
(2003). Putting the geometry all together leads to a ridge which 
is -0.45 pc long along the line of sight, -0.05 pc thick (width at 
half maximum of the density profile) and -0.42 pc long in the 
sky plane. Figure 9 compares the observed and modelled bright- 
ness profiles. The spatial shift between the maximum intensity 
positions observed at different wavelengths is well reproduced. 
The absolute brightnesses in the Herschel bands are also quite 
well reproduced by the model, except the 160 and 250 yum emis- 
sions which are underestimated. The modelled profiles are also 
narrower than the dat and, at shorter wavelengths, the PAH emis- 
sion is significantly overestimated by the model. 

7. Discussion 

Using our plane parallel model of the Orion Bar, assuming dust 
properties typical for the diffuse atomic ISM and with values 
of the local and column densities compatible with previous 
observations, we are able to reproduce very well the stratifi- 
cation of the dust emission inside the PDR seen in Herschel 
and Spitzer/IRAC data (Figs. 1 and 6). This indicates that the 
increase of the distance between the observed Bar and the 
illuminating star with increasing wavelength is a pure radiative 
transfer effect. This model is also able to reproduce at first order 
the global shape of the brightness profiles and the absolute 
brightness observed with Herschel. We conclude that a simple 
plane parallel model seems to be a reasonable approximation of 
the Bar. 

The discrepancy between the synthetic and observed brightness 
profiles in Herschel bands can be explained by two assump- 
tions we make. First, the PDR is modelled as a single slab 
without matter before and behind, and so the model does not 
reproduce the observed brightness outside the ridge. Second, 
the geometry may be not purely edge-on (e.g., Tielens et al. 
1993; Wen & O'DeU 1995; PeUegrini et al. 2009). Inclination 
effects, which are not taken into account in our study, can 
broaden the brightness profiles and the SEDs. Sub- structures 
which are ignored in our modelling can also broaden the profiles. 
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Fig. 9. Brightness profiles from Spitzer and Herschel observations (black), The photometric uncertainties are respectively 5% for 
IRAC, 20% for PACS and 7% for SPIRE (grey). The modelling using diff'use ISM dust properties (Tab 2) is shown by the red lines. 
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Fig. 10. Same as Figure 9 with a PAH abundance 7 times lower than in the diff'use ISM (Tab 2). The increase of brightness in the 
model at long wavelengths is due to a weaker absorption in the UV/visible. 
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Fig. 11. Same as Figure 9 with Ipdr = 0.25 pc. PAH abundance is a factor 3.3 lower than the canonical value in the diffuse ISM. The 
BG emissivity is enhanced by a factor 2 (Tab. 2). 

Table 2. Properties of the different models tested. 





7pAH {MIMh) 


^SamC 


iLamC 


^aSil 


^FIR 


/PDR(PC) 


Diffuse ISM model (Fig. 9) 


7.8 X 10-4 


1.65 X 10-4 


1.45 X 10-3 


7.8 X 10-3 




0.45 


PAH depleted model (Fig. 10) 


1.1 X 10-4 


1.65 X 10-4 


1.45 X 10-3 


7.8 X 10-3 




0.45 


PAH depleted model + 6bg enhancement (Fig. 11) 


2.36 X 10-4 


1.65 X 10-4 


1.45 X 10-3 


7.8 X 10-3 


2X6pV(^)^ 


0.25 



Notes. ^fir(^) emissivity presented in Fig.Al from Compiegne et al. (2011). 



The overestimated PAH emission in the Orion Bar while us- 
ing diffuse ISM dust abundances is not a surprise and could be 
due to several phenomena. 

1. Dust grains in the Orion bar undergo a very hard UV ra- 
diation field which can photo-destroy PAHs at the ioniza- 
tion front (Giard et al. 1994). This could be coupled with 
PAH sticking on BGs in the molecular cloud (Kim et al. 
1994). Reducing the PAH abundance by a factor 7, with- 
out changing any other parameter, allows a good match of 
the 3.6/im brightness profile (Fig. 10). Less PAHs reduces 
the UV absorption and therefore the emission in Herschel 
bands slightly increases (maximum brightness multiplied by 
factors 1.1-1.2 except at 160yum where it is unchanged). 

2. In a dense medium, where the radiation does not penetrate, 
dust grains are suspected to be coagulated (Tielens 1989), 
which increases the overall FIR emissivity epm by factors 
about 2-3 (e.g., Bernard et al. 1999; Stepnik et al. 2003; 
Paradis et al. 2009; Kohler et al. 2011). To be in agree- 
ment with the observations, a model with an increase of epiR 
must be counterbalanced by a decrease of the column density 



(through the Ipdr parameter) by the same factor. Taking into 
account this enhancement allows a good match of the bright- 
ness profiles, as illustrated in Figure 1 1 and a decrease of the 
length along the line of sight by factor 2-3 is still compatible 
with observational constraints. 
3. The length along the line of sight may be not constant. As 
mentioned above, in Herschel bands we probe the densest 
part of the Bar whereas at 3.6 yum we are sensitive to the sur- 
face emission which is very structured (first panel of Fig. 1). 
What appears as a bar in the sky plane could be a convex 
surface pointing to the illuminating star, and thus presenting 
a bigger depth along the line of sight in the densest part than 
near the ionization front. 

PDR models of the gas emission in the Orion Bar focus on 
the HII region and the zone near the ionization front, and 
fit the observations considering a length along the line of 
sight close to 0.15 pc (Pellegrini et al. 2009) while we take 
IpDR = 0.45 pc in our modelling. Therefore, the excess of 
the 3.6 yum emission in our model could be simply explained 
by a decrease of the length along the line of sight near the 
ionization front. 
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Fig. 12. Solid line: SED extracted from the radiative transfer 
model using diffuse ISM properties (Tab. 2) at the 70 yum peak 
position (Fig 9). Square: Integrated fluxes in Herschel bands. 
Dotted line: Modified blackbody fitted to the data. 

Einally, the observed absolute brightnesses at 160yum and at 
250yum are underestimated by factors 1.5 and 1.2, respectively. 
This comes from the dust model itself. With DustEM, optical 
properties of the BG component, made of silicates and large 
amorphous carbons, create an emerging spectrum with a value 
of the spectral emissivity index close to 1.55. Cooling and 
radiative transfer cannot change this. As a consequence, we fail 
to reproduce the SED of BGs, whatever the position on the Bar; 
grains are systematically cooler and spectra are systematically 
steeper in the data than in the model, as illustrated at the 70 [im 
peak position shown in Figure 12. SEDs of BGs at diff'erent 
position across the Bar (studied Sect. 4 and shown Fig. 5) reveal 
P values around 1.8 in the Bar. This discrepancy between the 
model and the data induces the underestimate of p seen in 
Figure 12. 



8. Conclusion 

We have presented the first observations of the Orion Bar pho- 
todissociation region using the Herschel Space Observatory in 
five broad bands (70 and 160 yum with PACS, 250, 350 and 500 
yum with SPIRE). The wavelength range and the spatial resolu- 
tion of Herschel allow us to map and resolve the emission of the 
biggest grains located in the densest zones of the Orion Bar. 
To extract SEDs, we have developed an original and efficient 
method to compute convolution kernels allowing to convert all 
maps to the same beam. Study of BG SEDs at diff'erent positions 
in the Orion Bar region reveals a cooling of dust (from 80 K to 40 
K) and an increase of the spectral emissivity index p (from 1 . 1 to 
2) as we go deep into the PDR. Combining the new Herschel ob- 
servations with ancillary data from Spitzer/IRAC gives us the op- 
portunity to probe the morphology of the Orion Bar from 3.6 yum 
to 500 yum. The global structure is the same at all wavelengths but 
we observe an increase in the distance star-Bar with increasing 
wavelength (shift of ~ 30" between the 3.6yum and the 500 yum 
profiles), which implies a stratification of dust properties within 
the Bar. 

We have used a radiative transfer model using typical dif- 
fuse ISM dust properties and abundance to model the brightness 
profiles observed with Herschel and Spitzer. This model repro- 



duces the stratification within the Bar using values of the column 
density and local density parameters that are in agreement with 
the observational constraints. This shows that this stratification 
is purely a radiative transfer eff'ect. 

However, using diff'use ISM dust properties and abundances 
in the model leads to several discrepancies in shape and absolute 
brightness between the modelled brightness profiles and the ob- 
servations. The assumptions we make about the geometry (con- 
sidering the PDR as a single slab, adoption of an edge-on geom- 
etry, without inclination and sub-structure eff'ects) can explain 
the diff'erences of width between the synthetic and observed pro- 
files. 

The relatively low absolute brightness at 3.6 yum might be ex- 
plained by PAH photo-destruction by the strong UV radiation, 
or by sticking of PAHs to the BGs in the pre-existing molecular 
cloud. A PAH abundance decrease by a factor 7 from the canoni- 
cal value in the diff'use ISM allows a good fit of the IRAC 3.6 jjm 
emission. A combination of a decreased PAH abundance (by a 
factor less than 7) with a decreased length along the line of sight 
is possible too. If so, an increase of the brightness in Herschel 
bands is needed to fit the data. Coagulation can explain such a 
behaviour since it can enhance the FIR emissivity by factors 2-3. 
An alternative is that the length along the line of sight is shorter 
at 3.6 yum since we probe the emission from the surface of dense 
structures. 
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